Thermal noise and dephasing due to electron interactions in non-trivial geometries 
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We study Johnson-Nyquist noise in macroscopically inhomogeneous disordered metals and give 
a microscopic derivation of the correlation function of the scalar electric potentials in real space. 
Starting from the interacting Hamiltonian for electrons in a metal and the random phase approxi- 
mation, we find a relation between the correlation function of the electric potentials and the density 
fluctuations which is valid for arbitrary geometry and dimensionality. We show that the potential 
fluctuations are proportional to the solution of the diffusion equation, taken at zero frequency. As 
an example, we consider networks of quasi-lD disordered wires and give an explicit expression for 
the correlation function in a ring attached via arms to absorbing leads. We use this result in order 
to develop a theory of dephasing by electronic noise in multiply-connected systems. 



I. INTRODUCTION 

Electronic noise generated by the thermal excitation of 
charge carriers has been observed and explained by John- 
son and Nyquist more than 80 years ago^and discussed in 
great detail in the literature since then. More recently, it 
has been found that this so-called Johnson-Nyquist noise 
is the main source of dephasing in mesoscopic systems 
at low temperatures of a few Kelvins, where phonons are 
frozen out. Dephasing puts an IR cut-ofF for interference 
phenomena, such as quantum corrections to the classical 
conductivity.'^l 

The current interest in this topic arises from stud- 
ies of dephasing in mesoscopic systems which consist of 
connected quasi-lD disordered wires, Fig. [T] including 
connected rings and grids.l^ It has been found (both 
experimentally^ and theoreticallj^sHIol) that dephasing de- 
pends not only on the dimensionality, but also on the 
geometry of the system. The noise correlation function 
is well-understood for macroscopically homogeneous sys- 
tems such as infinite wires or isolated rings, but has so- 
far not been studied in multiply-connected networks with 
leads attached at arbitrary points. The goal of this pa- 
per is to give a transparent and systematic description 
of the thermal noise properties for such systems. In par- 
ticular, we will derive an expression for the fluctuations 
of the scalar electric potentials for arbitrary geometries. 



Eq. (31), and a general expression for the correspond- 
ing dephasing rate, Eq. (45). Throughout, we assume 
that a description of the noise in terms of scalar poten- 
tials is sufficient, i.e. we neglect the fluctuations of the 
transverse component of the electromagnetic field (for a 
detailed discussion of the latter, see Ref. j2j). 

Let us start by reviewing simplified arguments to de- 
rive the noise correlation function: Johnson and Nyquist 
concluded that thermal noise in electrical conductors 
is approximately white, meaning that the power spec- 
tral density is nearly constant throughout the whole fre- 
quency spectrum. If in addition the fluctuations are un- 
correlated for different points in space, a correlation func- 
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FIG. 1. A network of wires: We are interested in the noise 
correlations between arbitrary points x and y of multiply- 
connected networks attached to leads (denoted by the usual 
ground symbol) at arbitrary points. 



tion for the random thermal currents in the classical limit 
is independent of frequency uj and momentum q. The 
power spectrum of the current density reads 



■}j\^){q,Oj) = 2Tao 



(1) 



Here, ctq — e^vD is the Drude conductivity of the disor- 
dered system, D and v are the diffusion constant and the 
density of states, respectively. Naively applying Ohm's 
law, j{q) = (ToE{q), to Eq. ([T]) and using the rela- 
tion between the electric field and the scalar potential, 
eE{q) — ~iqV{q), we find 



{\V\')iq,u:) 



2Te^ 1 
CTo 



(2) 



The correlation function, Eq. ([2]) , corresponds to the cou- 
pling of agiven electron to the bath of the surrounding 
electrons .El Thus (I^P) describes the process of succes- 
sive emission and re-absorption of a photon, which is de- 
scribed effectively by the scalar potential V. The factor 
l/<7^ coincides with the solution of a diffusion equation 
in an infinite system, which reflects the fact that the cur- 
rents, Eq. are uncorrelated in space. 

These simple arguments are based on the homogeneity 
of the system and have assumed a local relation between 
potential and current, whereas transport properties in 
disordered metals are substantially non-local. In this 
paper, we derive an analogy of Eq. ([2| for disordered sys- 
tems with arbitrary geometry and dimensionality; this 
will in particular apply to networks of disordered wires. 
A detailed calculation, which takes into account all prop- 
erties of the mesoscopic samples, has to be done in the 
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real-space representation. Starting points are the usual 
linear response formalism and the fluctuation-dissipation 
theorem (FDT) ^ Although most ingredients of the fol- 
lowing discussion will be familiar to experts, we hope that 
the manner in which they have been assembled here will 
be found not only to be pedagogically useful, but also 
helpful for further theoretical studies. 

The paper is organized as follows: In Section|llj we pro- 
pose a heuristic description of the potential fluctuations. 
In Section |III[ we review a microscopic approach to the 
noise correlation function, based on a relation of the fluc- 
tuations of the scalar potentials to the fluctuations of the 
density, using the random phase approximation (RPA). 
In Section |IV] we evaluate the density response function 
X for disordered systems by using a real-space represen- 
tation for arbitrary geometries. We apply this result to 
the noise correlation function in Section |Vl In Section 
|VI[ we show how the noise correlation function can be 
calculated for networks of disordered wires. Finally, in 
Section IVIII we discuss the relation to the fundamental 
problem of dephasing by electronic interactions. 



II. HEURISTIC DESCRIPTION OF 
POTENTIAL FLUCTUATIONS 

A description of fluctuations in metals within the linear 
response formalism naturally starts with an analysis of 
the density fluctuations in the model of non-interacting 
electrons described by the standard frcc-clcctron Hamil- 
tonian This system is perturbed by an external 

scalar potential V{x,t) coupled to the density operator 
h{x): 

= J dxV{x,t)n{x) . (3) 
The response of the (induced) charge density. 



where 



dte" 



= - Jdyx{x,y,uj)V{y,uj) , (4) 
is governed by the (retarded) density response function: 

POO 

x{x,y,uo) = I / dte'(-+»")* ([n(a;, <), n(y, 0)]) . (5) 



Here (• • ■)pert ^'^'^ (' ' ') denote quantum/statistical av- 
eraging with respect to the perturbed and unperturbed 
Hamiltonian, respectively. The FDT relates the equilib- 
rium density fluctuations to the imaginary (dissipative) 
part of the response function, 

/oo 
dte^"'(n(a;,i)n(y,0)) (6) 
-OO 

= F(a;)Ini[x(a;,y,w)] , (7) 



1 - e-'^IT 



(8) 



In writing Eqs. (|6])-(|8]) we have exploited detailed bal- 
ance and time-reversal symmetry. The latter implies 
X{x,y,uj) = x{y,x,uj)^ 

The question which we are going to address in this pa- 
per is how to characterize the fluctuations of the electric 
potential V . For this purpose we consider the "dual" 
case, where some external density nc-x_t{x,t) is the per- 
turbation that couples to the "potential operator" 



^ JdxV{x)n,^t{x,t). 



(9) 



The linear response of V to the perturbation can be writ- 
ten as 



(F(a;,a;))port J dy T(a;, y, cj) ncxt(2/, w) , (10) 

defining the response function T. In analogy to Eq. ([7|, 
the response function also characterizes the equilibrium 
fluctuations of the potentialfl^ 

{\V\^)ix,y,Lu) - F{oj)lm[-T{x,y,L,)] . (11) 

Calculating the response function T{x,y,uj) is a com- 
plicated task because it requires precise knowledge of the 
potential operator V{x). Instead, we can identify the po- 
tential V(x,uj) in Eq. Q with the response {V{x,uj))pc,-t 
in Eq. ( 10 ) to relate T to x- In the limit of strong 



screening in good conductors (called the unitary limit^), 
electroneutrality is satisfied locally. Therefore, the in- 
duced charge exactly compensates the external charge: 
nind(a',w) = — next(a3,w). Now inserting Eq. Q into 
Eq. (10) (or vice versa), we obtain 



J dx' T{x, x', uj)xix', y, w) = S{x - y) . 



(12) 



If X is known, Eqs. (11) and (12) allow one to calculate 



the correlation function of the scalar potential. 

Let us recall the well-known case of macroscopically 
homogeneous diffusive systems. The express ion fo r the 
disordered averaged response function x rea.dff^^^ 



xiq,uj) = V 



Dq^ — iuj 



^ = 1/T(q,w), 



(13) 



where we used Eq. ( 12 ). Inserting Eq. ( 13 1 into Eq. (Ill, 
we find 



{\V\'){q,oj) = 



1 ojF{u) 



V Dq^ ' 

which reduces to Eq. ^ in the limit uj -^T. 



(14) 



3 



UKPA{x,y,uj) = X 

= aS'v/x/x/w ?7 + X 

Uo 

FIG. 2. The Coulomb interaction in the RPA according to 
Eq. ITfl. 
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III. NOISE CORRELATION FUNCTION FOR 
ARBITRARY GEOMETRIES: MICROSCOPIC 
APPROACH 

In Eqs. ([3]) and ([9|, we introduced the operators n and 
V assuming that either V{x,t) or next{x,t) are external 
perturbations. In fact, the fluctuations originate inside 
of the system and the starting point of a microscopic de- 
scription is the part of the Hamiltonian, which describes 
electron interactions, 



da;dy Uq{x, y)-ip\x)-i}:\y)'ip{y)i)[x) , (15) 



where Uo{x,y) is the bare Coulomb interaction. In the 
mean-field approximation, Eq. ( 15 ) gives rise to a correc- 



tion, called Hartree contribution, to the electron energy: 

^^(Hartrec) ^ / dxdy Uo{x , y) {h{x)} {h{y)} , (16) 



where h{x) = {x)i^{x). 

The Coulomb interactions are dynamically screened, 
which can be accounted for in the framework of the RPA, 
provided that the electron density is high, 

URPAix, y,u}) = Uoix,y) - J dx'dy' Uo{x,x')x{x' ,y' ,u}) 

xURPA{y\y,uj), (17) 

see Fig. [2j Note that — x [see the definition in Eq. ( [s])] i s 
equal to the bubble diagrams of Fig. [2] see e.g. Ref. |l8j. 
In Appendix |Aj we recall how to obtain Eq. (171 within 
a self consistent treatment of the screening problem. 

Using the RPA in Eq. ( 16 ) and comparing the result 
with equation Eq. (|9| , we observe that the potential fiuc- 
tuations are due to electronic interactions and that the 
operator of the scalar potential is given by 



(18) 



y(a;,w)= dyUviPAix,y,u})n{y) 



Eq. ( 18 1 allows us to relate the correlation function of 



the potentials to the correlation function of the density 
fiuctuations: 

{\V\^)ix,y,Lj) = Jdx'dy'URPAix,x',Lj) 

X (|n|2)(a;',y',c^)C/i^pA(y',y,c^). (19) 

By inserting Eqs. (17) and ([t]) into Eq. (19), re-ordering 
the terms in the RPA series and using the fact that Uq 
is real, we find (see Fig. [3]) 

(|y H {x, y, u;) = Fiu) Im[-C/RPA(a;, y, u)] . (20) 




= Im[-I7RPA 



FIG. 3. Diagrammatic proof of Eq. ( |20[ ) by using Eq. ([19| 
and Eq. ^ (i.e. Fig.[2|. 



We emphasize that the derivation of Eq. ( 20 ) has not used 



any other assumption than the RPA. Thus, Eqs. ( 17 1 and 



(20) are a microscopic (and more rigorous) counterpart 



of the phenomenological Eqs. (11) and (12) 



IV. DENSITY RESPONSE IN DISORDERED 
SYSTEMS: CALCULATIONS IN COORDINATE 
REPRESENTATION 

In disordered metals the motion of the electrons is dif- 
fusive, provided that kp^ ^ i L where kp is the 
Fermi wave-vector, £ the mean free path and L the sys- 
tem size. It can be accounted for by substituting the 
disorder-averaged density response function, x, into the 
phenomenological Eqs. (11) and (12) or the microscopic 
Eqs. (17) and (20). The function x has been calculated 
for macroscopical y homogeneous systems by Vollhardt 
and Wolfle.-^- In the following, we will show how to gen- 
eralize their calculation to inhomogeneous systems. A 
useful starting point is a coordinate representation of the 
density response function, Eq. ([5]), in terms of the ad- 
vanced and retarded Green's functions, G^/'^{x,y,uj), 

X{x, y,u)^^Jde[ [fie + lo) - fie)] (21) 
X G«(a?,y,e + w)G^(y,a;,e) 



+ fie)G'^ix,y,e + uj)GRiy,x,e) 

- fie + Lo)G^ix,y,e + uj)G^iy,x,e)'^ , 

see Fig. |4](a). Here, /(w) is the Fermi distribution func- 
tion, and denotes disorder averaging. The combi- 
nations G^G^ and G^G^ give short-range contribu- 
tions, since the average of the products decouple, e.g. 
G«G« ~G^-G^ + Oil/kpi), and the disorder aver- 
aged Green's functions G^ and G^ decay on the scale 
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X XX 

T{x,y,uj)^ = 6{x - y) + ^^y 

X x x' 



FIG. 4. (a) Equation for the disorder averaged density re- 
sponse function; solid lines denote the disorder averaged re- 
tarded and/or advanced Green functions, cf. Eq. (21 1, (b) 
Equation for the impurity vertex; the dashed line represents 
impurity scattering, cf. Eq. (23 I. 



£ <^ L. We will consider contributions to the thermal 
noise which are governed by distances larger than £, cf. 
Ref. [2] ■ Therefore, details of the behavior on short scales 
are not important for our purposes and we replace the 
short-range contributions, G^G^ and G^G^, by a delta 
function. The long range contributions G^G^ can be 
calculated by standard methodSjIi^ 



G^{x, y, e + uj)G^{y, x,e) ^ J dx' G^{x, x\ e + w) 

y.G^{x\x,e)T[x',y,u), (22) 

where r(a;,y,cj) is the impurity vertex function, 

1 



r(a;,y,w) = 5{x-y 



27ri/T 

X G^{x\x,e)T{x' ,y,uj 



Ax' G'^ix,x',e + uj) 



(23) 



(the factor l/27rz/T, where r = i/vp is the transport time, 
originates from the impurity line), see Fig. Qb). The 
short-ranged product, G^ ■ G'^, can be expanded as G^ ■ 
qA ^ 2tti't S{x — x') [l+iujT + TDAx\ , which is obtained 
by transforming the product to momentum space and 
expanding in the transferred momentum q and frequency 
w, realizing that terms of order q vanish due to symmetry. 
As a result, Eq. (23 1 reduces to a diffusion equation: 



{-iuj-DA^)T{x,y,uj)^-S{x-y), (24) 



where D = vp^/d is the diffusion constant for a d dimen- 
sional system. Thus, the vertex function is proportional 
to the diffusion propagator, r{x,y,Lu) = P{x,y,Ld)/T. 

Collecting the short- and long-range contributions and 
taking the limit T ^ e^?, we obtain from Eq. (21 ) 



X{x, y,uj) = v {5{x -y) + iujP{x, y, uj)) 



(25) 



Eq. ( 25 1 is valid for arbitrary geometries since it is based 



only on the diffusive approximation and does not require 
macroscopic homogeneity. 



V. NOISE CORRELATION FUNCTION IN 
DISORDERED SYSTEMS 



Let us simplify Eq. (17) for a disordered conductor. 
Using Eq. (l25|) and 



1 

47re2 



AxUq{x, x') = 5{x - x') , 



(26) 



Eq. (17) can be written as 



1 — U^pp^{x, y, Lo) j Ax' P{x, x' , lo) 

X ;7RPA(a;', y, uj) = -6{x - y) , 



(27) 



where we intr oduced the Thomas-Fermi screening wave- 
vector K = V 47re^z/, which corresponds to the inverse 
screening length in three dimensional (3D) bulk systems. 
The kernel of Eq. (271 is a solution to the diffusion equa- 



tion ( 24 ) which can be expanded in terms of eigcnfunc- 



tions of the Laplace operator. Consequently, the kernel 
is always separable and Eq. ( |27| has a unique solution 
(see e.g. Ref. [12] for details on how the solution can be 
found). Using the semi-group property of the diffusion 
propagators. 



dx' P{x,x',uj)P{x',y,0) 
one can check that 



[P{x,y,0) - P{x,y,u;)] , 



UnPAix,y,uj) = 



1 1 



P{x,y,{)) 



(28) 

satisfies Eq. (27 1. In practice, the 3D Thomas- Fermi 
screening length is a microscopic scale, thus the 



typical value of the first term of the rhs. of Eq. (28 1, 



(ZJq^yp — iwtyp)^^, is larger than — l/AiraQ for 

good conductors (this is the so-called unitary limit, for 
details see Ref. [9]): 



> 



I^Qtyp - ^'^typl Dk 



(29) 



In this limit, using the diffusion Eq. (24), we obtain from 
Eq. (l28l): 



1 



Ukpa{x, y,uj)^ - [5{x - y) - iujP{x, y, 0)] . (30) 

We remind that P{x,y,0) is always real. As a result, 
Eqs. (pO| and (pOl) yield 



(|U|2) [x, y,uj)^^ ujF{u:)P{x, y, 0) , (31) 

where F{uj) is given by E q. ([s] ). The real-space demon- 
stration of Eqs. (30) and (13 ip for macroscopically inho- 



mogeneous systems, are among the main results of the 



paper. It is worth emphasizing the frequency-space fac- 
torization of the correlator, which plays an important role 
in the theory of dephasing, cf . Section |VII| The relation 
of Eq. ( [31] ) to the correlation function of the currents, 
Eq. ([1]), is discussed in Appendix [b| and allows to put 
the presentation of the introduction on firm ground. 

Note that Eq. (29) allows one to neglect the term 
Is.x/k^ in Eq. ([27| and, thus, to reduce Eq. ([27|) to the 



form of the phenomenological integral equation ( 12 1, with 
C/rpa taking the place of T. [The same replacement leads 
from Eq. (11) to Eq. (20).] In other words, the electric 



potential of the fluctuating charge densities itself is neg- 
ligible when screening is strong enough (i.e. good con- 
ductors in the unitary limit), justifying a posteriori our 
assumptions in the phenomenological Section [Til 

The fact that the correlation function of the potential 
is proportional to the solution of the diffusion equation 
at zero frequency, cf. Eq. (31), may be understood as a 



relation function, Eq. (31), for a network of disordered 
wires. The main ingredient to Eq. (31) is the solution 



of the diffusion equation (24) at zero frequency. Wires 



allow a quasi- ID description of diffusion, where trans- 
verse directions can be integrated out since P{x,y,ui) is 
assumed to be constant on the scale of the width of the 
wire. As a result, we replace P{x,y,uj) — > P{x,y,uj)/s, 
where s is the cross-section of the wires and P{x, y, uj) 
solves the ID diffusion equation in the network, x and 
y being coordinates along the wires. Recently, effective 
methods have been developed to solve the resulting dif- 
fusion equation for arbitrary networks. ^"^^^ We will 
review these methods in this Section and evaluate the 
noise correlation function for a simple example. 

We start by introducing some basic notations: A net- 
work is a set of vertices, labeled by an index a, connected 
via wires of arbitrary length, say Zq,^ for the wire connect- 
ing vertices a and B. Let us define a vertex matrix A4 




non-local version of the Johnson-Nyquist theorem, since 
P{x,y,0) can be related to the classical dc-resistance 
TZ{x,y) between the points x and y (see Ref. j20j): 

n{x, y) - ^ I ^ [P{x, X, 0) + P{y, y, 0)] - Pix, y, 0)| . 

For example, in an infinitely long quasi-lD wire of 
cross-section s, the solution of the diffusion equation 
is P{x,y,0) — —\x — y\/{Ds) where x is the compo- 
nent of X along the wire. Hence we recover a resis- 
tance proportional to the distance between the points, 
n{x,y) ^\x- y\/isao). 



VI. NOISE CORRELATION FUNCTION IN 
NETWORKS OF DISORDERED WIRES 

Let us now illustrate the calculation of the noise cor- 



FIG. 5. The network corresponding to a symmetric ring made 
of four wires connected to two absorbing leads. The length 
of the arcs is / and the length of the connecting arms is b. 
Vertices are labeled by numbers a — 1,2,3,4. Vertices "1" 
and "4" denote the points where the ring is connected to the 
arms. Vertex "2" is always placed in the upper arc, defin- 
ing the running coordinate x. Vertex "3" determines the y- 
coordinate and is placed either in the adjacent arc (panel a) 
or in the left arc (panel b) or in the same arc (panel c). 



where a^^ = 1 if the vertices a and /3 are connected and 
Oa/? = otherwise. The solution of the diffusion equation 
at zero frequency between arbitrary vertices a, j3 of the 
network is given by the entries of the inverse matrix Ai 
divided by the diffusion constant: 



P(a,/3,0) - {M-^)o.p/D. 



(34) 



This allows us to calculate the noise correlation function 
between arbitrary points of a network by inserting ver- 
tices and inverting M. As an aside, note that arbitrary 
boundary conditions can be included in this scheme eas- 
ily (see Refs. i20]and[22] for details). 

Let us consider the network shown in Fig. [5] represent- 
ing a ring connected to absorbing leads. For simplicity, 
we assumed that the ring is symmetric: the two arcs are 
of the same length I and the connecting arms of length 
b. We evaluate the noise correlation function for two 
points in this network by inserting two vertices, called 
"2" and "3". Vertex "2" is always placed in the up- 
per arc, encoding the running coordinate x in the length 
of the connected wires. Vertex "3" determines the y- 
coordinate and is placed either in the lower arc or in the 
left connecting arm or in the upper arc. In the first case. 
Fig. [5ja), the vertex matrix, Eq. (33), is given by 



/ 







\ 
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FIG. 6. (color online) The solution to the difTusion equation 
at zero frequency, P{x, y, 0) oc {\V\^){x,y,uj), where (|V^|^) is 
given by Eqs. (36l-(38l, for a fixed coordinate x in the upper 



arm of the ring (indicated by the dot), as a function of y 
traversing the network. {\V\^){x,y,uj) is linear in y and its 
derivative has a discontinuity at y = x. 



At low temperatures, dephasing is dominated by elec- 
tron interactions, that can be accounted for through 
a contribution to the phase accumulated by two time- 
reversed interfering trajectories in a fluctuating electric 
fielcP. 



$[a;(T)] = [ dr [V{x{t),t) ~ V{x{T),t - t)] 
Jo 



(39) 



When averaged over the Gaussian fluctuations of the elec- 
tric field, (e**)y — e~2(* )v yields a phase difference 
which cuts off the contributions of long electronic trajec- 
tories. Introducing the trajectory-dependent dephasing 
rate r[x{T)] = ^ ($[a;(r)]^)y, the weak locahzation cor- 
rection takes the forml^MSHHl 



The diffusion propagator is then given by P{x, y, 0) = 
{M.~^)2?, I D , and we obtain from Eq. (31 ) the correlation 
function as a function of the running coordinates x, y 
e [0,?]: 



{\V\''){x,y,u) = 



ujFiuj) b{l{2b + l)-{x + y)l + 2xy) 



Dvs 



l{4b + l) 



(36) 



When vertex "3" is placed in the connecting arm, x G 
[0,/] and 2/ G [0,6] [Fig. [sjb)], we get 

(37) 



Di^s 



Finally, when vertex "3" is placed in the same arc of the 
ring as vertex "2" [see Fig.[5jc)], following the same logic 
we obtain, with < x < y < I, 

{\V\')ix.y,u;) 

ujF{uj) bl{2b + l)+ xl{3b + 1)- ybl - xy{2b + I) 



Dvs 



l{Ab + l) 



(38) 



All other configurations can be found by symmetry ar- 
guments. We plot P(x,y,0) for y traversing the whole 
network in Fig. [6] Note that the resulting function is 
linear in y and its derivative has a discontinuity at y = a; 
(cf. Ref. [2]). 



VII. APPLICATION TO DEPHASING 

The precise characterization of potential fluctuations 
is very important in studying phase coherent properties 
of disordered metals at low temperatures. To be specific, 
let us discuss a particular coherent property: the weak 
localization correction to the conductivity. Let us recall 
that the weak localization (WL) correction Act = ct— ctq is 
a small contribution to the averaged conductivity arising 
from quantum interference of reversed diffusive electronic 
trajectories.!^ 



Aa{x) 



2e^D 



dt P{x, X, t) (e 



-tr[x(r)]. 



{x(t)} 



where 



(40) 

){x(t)} is the average with respect to closed 
diffusive trajectories of duration t starting from x (not 
to be confused with the thermal average (• • ■)y over the 
electric potential V). The phase fluctuations can then be 
related to the potential fluctuations: 



^-mxf)= J^drdr' 



duj 
. 2^ 



x{\V\'}^{x{r),x{r'),uj). 



Here we have introduced a new noise correlator, 



{\V\'')^{x, y,io) = - ujF^{uj)P{x, y, 0) , 



(41) 



(42) 



obtained from Eq. ( 31 1 by replacing F{uj) with a modified 
function Fip{uj) (given below), on the origin of which we 
now comment. Equation (20) is well-known in the the- 



ory of dephasing: its version symmetrized with respect to 
frequency arises naturally when compa ring t he diagram- 
matic calculation of the dephasing time^SEll with the in- 
fluence functional approach d escribin g electrons moving 
in a random Gaussian field y IMMin] Diagrammatically, 
the symmetrized Eq. (20| represents the Keldysh com- 
ponent of the screened electron interaction propagator, 
the only substantial difference being that the diagram- 
matically calculated correlation function involved in the 
dephasing process acquires so-called "Pauli-factors" that 
account for the fact that the Fermi sea limits the phase 
space available for inelastic transitions.!^ These factors 
lead to the following replacement of the function F{uj) in 
Eq. (IM and here also in Eq. (fsTl); 



Fie.] 



-.coth(./2T)^^.^^,^^^^^^.^^,(.) 



(43) 

This restricts the energy transfer to |a;| < T^^^hni does 
not affect the factorization of the correlator. Inserting 
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Eq. Q into Eq. ^ leads to 
1 .0. 2T 



drdr' P{x{t),x{t'),0) 



-oo 27r 



g-icj(T-T') _ g-ja;(T+r'-t) 



2T 



^^^H- (44) 



The fact that the frequency dependent function ■^Fip{uj) 
is symmetric aUows us to add to P{x{t),x{t'),0) the 
term -1[P{x{t),x{t),0) + P{x{t'),x{t'),0)], which 
does not contribute to the integral (44). Therefore we 



finally end up with the following expression for the de- 
phasing rate, 

T[x{T)]^e^T f%- fdr' [Srir + r' - t) - S^r - r')] 
Jo ^ Jo 

xn{x{T),x{T')) , (45) 



written in terms of the resistance TZ, defined in Eq. ( 32 1 



The function drit), a broadened delta function of width 
1/T and height T, is the Fourier transform of ^F,^{uj), 
which is given by 



SriT) = ttT w{ttTt) 



ycothy - 1 

^(y) = r-r2 ■ (46) 

smn y 



Eq. (45), which is one of the main results of our paper, 
generalizes the results obtained in Refs. [5J [51 [53J and 
128] for an infinite wire and an isolated ring to arbitrary 
geometry. In the classical noise limit, T — > oo, 6t{t) 
may be replaced by a (5(r) function: The second term of 
Eq. (45) vanishes and we recover the results of Refs. ^ 
andllOr 



Let us now illustrate Eq. ( 45 1 by calculating the de- 



phasing time for the well-understood case of one and 
two-dimensional isolated simply-connected samples. The 
dephasing time can be extracted from the condition 



1 ^ r(r^ 



(47) 



where r(t) is given by the functional Eq. (45), averaged 
over the typical closed random walks x{t) of duration t 
in the system. The problem is governed by the interplay 
of three time-scales: The Thouless time rxh = L^/D, 
depending on the system size L, the thermal time tt = 
1/T (related to the thermal length Lt = \/D/T), and 
the dephasing time r^p. 

(i) Diffusive regime, tt <C r^^ <C rxh {Lt <C <C L): 
this is the regime considered in Refs. P and , where 
the width of the broadened delta functions in Eq. ( 45 ) , 



tt, is the shortest time-scale. Thus, when averaging over 
paths £c(t), the characteristic length-scale \x — y\ entering 
the resistance TZ{x,y) can be determined as follows: For 
the first 5t term this length is governed by free diffusion, 
since \x{t) — x{t — t)\ ^ VDt, hence jx — y| ~ V Dr. 
For the second term the characteristic length is set by the 
width of the delta function, \x—y\ ^ \JDtt- In ID, where 
TZix, y) ~ |x— j/I/ctqs, the first term dominates and we im- 
mediately obtain from Eq. ((47| l/r^ {e'^y/DT/aos)'^^^. 



In 2D, the diffuson at zero frequency is logarithmic as well 
as the resistance (32), TZ{x,y) ^ \n{\x — y\)/aod, where 
d is the width of the sample, which can be understood 
from the fact that the resistance of a plane connected at 
two corners scales logarithmically with the system size. 



Eq. (|47]) gives 1 e'^TT^\ii(TT^)/aod, and for the de- 
phasing time: ~ e^Tln.{e^/aod)/aod. 

(ii) Ergodic regime, tt <C rxh < (Lt < L < L^): 
the width of 6t (t) in Eq. ( 45 ) is still the shortest time- 



scale but, in contrast to (i), the typical trajectories x{t) 
explore the whole system, setting the length-scale of dif- 
fusion to the system size L, cf. Refs. and [7]. In 
full analogy to the diffusive regime, but replacing ^/Dt 
by L, we find for ID, 1/t^~ e'^LT/aos, and for 2D, 
l/r^ ^ e^Tln(rTh/'''T)/o'orf-'^° These examples show that 
for non-trivial geometries dephasing due to electron in- 
teractions cannot be accounted for through a unique de- 
phasing rate depending only on dimensionality, but must 
be described by a functional of the trajectories x{t) since 
the qualitative behavior of r^p follows from the geometry 
dependent typical distance \x{t) — x{t')\. 

For sufficiently low temperatures, on the other hand, 
Eq. (45) is capable to describe the crossover to a OD 



regime, where, apart from a dependence on the total sys- 
tem size, geometry becomes unimportant: 

(iii) OD regime, rxh < tt < r^^ {L ^ Lt <. L^): 
here, the width of the delta functions in Eq. (45), t^, 
is larger than Txh- Hence, the trajectories reach the er- 
godic limit x{t > Txh) ~ L before the electric poten- 
tial has significantly changed: Dephasing is strongly re- 
duced. Let us denote the maximal resistance reached at 
the ergodic limit as 7?.erg and replace the resistance in 



Eq. (45 ) by 7?, — 7?. — 7?.erg, without changing the result. 



since 7?.org is constant and its contribution vanishes af- 
ter integrating over r and r'. The difference TZ — 7?.crg 
is nonzero only during time differences t — t' < rxh, 
before reaching ergodicity. Thus, the leading contribu- 
tion comes from the second 5t term in Eq. ( |45[ ), which 
is constant at its maximum T during such short time- 
scales. We find 1 e^rVy /J'^^dr [7^(a;(T),0) - H. 

and since the TZ, 
ing time l/r^p 



I. J- 'yJO V'^\'^\' It^J '^crgj 

org term dominates, we obtain a dephas- 
"^T^TTh'T^-erg, independent of geometry 



and with the characteristic ^ behavior 



VIII. CONCLUSIONS 

In this paper we have considered fluctuations of the 
scalar electric potentials in macroscopically inhomoge- 
neous metals. We have shown how to relate the den- 
sity fluctuations to the potential fluctuations, emphasiz- 
ing the role of electronic interactions, provided a real 
space derivation of the density response function and il- 
lustrated these general ideas for the case of networks of 
metallic wires. Finally we have obtained a trajectory- 
dependent functional, Eq. (45), which describes dephas- 



ing by electron interactions for arbitrary geometries and 
accounts for the quantum noise contribution. When ap- 
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plied to networks, Eq. ( 45 ) can describe the full crossover Convolution with the Coulomb interaction gives Eq. (171 



from the OD to the ID and the 2D regime. 
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Appendix A: Self consistent analysis of screening 



We recall here how to obtain Eq. (17) using a self- 
consistent treatment of screening in real space.f^^l Starting 
points are the following three equations: (i) the excess 
charge density is decomposed into external and induced 
contributions 



Sn{x,uj) = ncxt(a;,w) + nind{x,uj) 



(Al) 



(ii) The induced charge is related to the potential V{x, u) 
by the density response function, cf. Eq. (|4]): 



Jdy x{x,y,oj)V{y,oj). (A2) 



(iii) The Poisson equation 

A V{x, uj) = -47re^ 6n{x, oj) . 



(A3) 



Self-consistency lies in the fact that the response involves 
the screened potential V{x,uj) and not the bare "exter- 
nal" potential related to next(a3,a;). The screened effec- 
tive interaction between electrons U^pa{x , y , lo) is ob- 
tained by placing an external charge at y, so that the 
external density is ncxtix,uj) — 5{x — y), and associating 
the resulting screened potential V{x,ijj) in Eq. (A3) with 
t^RPA(a;,y,a;). We obtain 



1 



47re 



2^xt^RPA(a',y,w) 



j^x' x{x,x',uj)URPA{x',y,uj) ^5{x-y). (A4) 



Appendix B: Current density correlations 

We discuss here the relation between the density and 
the current density correlations. The response of the (in- 
duced) current density is characterized by the conductiv- 
ity tensor cr, 

(ja(a^, w))ncq = j dy cTapix , y , uj) Ep{y , Uj) , (Bl) 
which is related to Eq. ([5| by current conservation: 



^a^pCTapix, x\ Uj) = -iuje'^x{x, x' , Oj) . 



(B2) 



The thermal fluctuations of the current density can be 
obtained from {jajl){x,y,uj) = wF(cj) Re[tTa^(a;, y, w)], 
in analogy to the discussion in Section assuming time- 
reversal symmetry, aapix^y^uS) = ai3a{y,x,uj). 

Let us now examine the case of disordered metals. 
The classical contribution to the averaged non-local dc- 
conductivity has been derived in Ref. [T^]. Their result 
can be generalized straightforwardly to non-zero frequen- 
cies, 

CTa/3(a;,a;',w) = CTo [5ai)5{x - x') - a^'pP{x,x' ,u;)\ , 

(B3) 



which obeys the condition (B2) with Eq. ([25|) substituted 
for X- For the current correlations, we find—- 



{jajl){x,x',uj) =aoujF{uj) {Sai3S{x - x') 

- D\/cVpRe[P{x,x',uj)] } 



(B4) 



Since the diffuson, P{x, x' , lS), decays exponentially on a 
length scale L^^ = ^JD/lo^ this expression shows that 
current correlations can be considered as purely local 
over the scale ||a; — x'\\ ^ L^, i.e (|jp) (a;, a;', oj) ~ 
aQijjF{ijj)6{x — x'). In the limit of classical noise, 
F{uj)u} ~ 2T, we recover precisely Eq. ([!]). 
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